Search for: All records

Abstract While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete matter from a surrounding disk. To study the dynamics of this system requires simultaneously describing the evolving spacetime and the magnetized plasma. We present the first relativistic calculation simulating two equal-mass, nonspinning black holes as they inspiral from a 20M(G=c= 1) initial separation almost to merger. Our results imply important observational consequences: for instance, the accretion rate $\dot{M}$onto the black holes first decreases and then plateaus, dropping by only a factor of ∼3 despite the rapid inspiral. An estimated bolometric light curve follows the same profile, suggesting some merging SMBBHs may be significantly luminous past the predicted circumbinary disk decoupling. The minidisks are nonstandard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, “sloshing” streams transfer mass between minidisks, carrying kinetic energy at a rate sometimes as high as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives time-varying minidisk tilts. These complex dynamics all contribute to unique cyclical behavior in the light curves of late-time inspiraling SMBBHs. The poloidal magnetic flux on the black holes is roughly constant at a dimensionless levelϕ∼ 2–3, but doubles just before merger; for significant black hole spin, this flux predicts powerful jets with variability driven by binary dynamics, another potentially unique electromagnetic signature. This simulation is the first to employ our multipatch infrastructure PatchworkMHD, decreasing the computational expense to ∼3% of conventional single-grid methods’ cost. 

Abstract Accretion of debris seems to be the natural mechanism to power the radiation emitted during a tidal disruption event (TDE), in which a supermassive black hole tears apart a star. However, this requires the prompt formation of a compact accretion disk. Here, using a fully relativistic global simulation for the long-term evolution of debris in a TDE with realistic initial conditions, we show that at most a tiny fraction of the bound mass enters such a disk on the timescale of observed flares. To “circularize” most of the bound mass entails an increase in the binding energy of that mass by a factor of ∼30; we find at most an order-unity change. Our simulation suggests it would take a timescale comparable to a few tens of the characteristic mass fallback time to dissipate enough energy for “circularization.” Instead, the bound debris forms an extended eccentric accretion flow with eccentricity ≃0.4–0.5 by ∼two fallback times. Although the energy dissipated in shocks in this large-scale flow is much smaller than the “circularization” energy, it matches the observed radiated energy very well. Nonetheless, the impact of shocks is not strong enough to unbind initially bound debris into an outflow. 

Abstract We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with a minidisk around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by a m = 1 overdensity called the lump. We compare our new spinning simulation with a previous non-spinning run, studying how spin influences the minidisk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the minidisks. We find that minidisks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases we find that most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes, and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.